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In this work a phenomenological Stochastic Differential Equation (SDE) is proposed 
for modelhng the time-evolution of the radius of a pre-critical molecular cluster dur- 
ing nucleation (the classical order parameter). Such a SDE constitutes the basis for 
the calculation of the (nucleation) induction time under the Kramers' theory of ther- 
mally activated escape processes. Considering the nucleation stage as a Poisson's 
rare-event, analytical expressions for the induction time statistics are deduced for 
both steady and unsteady conditions, the latter assuming the semiadiabatic limit. 
These expressions can be used to identify the underlying mechanism of molecular 
cluster formation (distinguishing between homogeneous or heterogeneous nucleation 
from the nucleation statistics is possible) as well as to predict induction times and in- 
duction time distributions. The predictions of this model are in good agreement with 
experimentally measured induction times at constant temperature but agreement is 
not so good for induction time statistics. Stochastic simulations truncated to the 
maximum waiting time of the experiments confirm that this fact is due to the time 
constraints imposed by experiments. Correcting for this effect, the experimental and 
predicted curves fit remarkably well. Thus, the proposed model seems to be a versa- 
tile tool to predict cluster size distributions, nucleation rates, (nucleation) induction 
time and induction time statistics for a wide range of conditions (e.g. time- dependent 
temperature, supersaturation, pH, etc.) where Classical Nucleation Theory (CNT) 
is of limited applicability. 
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I. INTRODUCTION 

The initial (nucleation) stage of a first-order phase transition, during which the precursors 
of the new phase appear as a result of random density fluctuations, still remains a challenging 
problem in the field of condensed matter. Throughout the last century, numerous theoretical 
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models have been developed with the aim of describing such phenomenon. These can be 
grouped in two categories depending on how the energy needed to form a molecular aggregate 
is derived: i) the cluster approach, originally introduced by the pioneers of the nucleation 
theorjff^ and ii} the density-functional approach (DFA), which was applied to nucleation 
first by Cahn and HilliarcP. 

In the cluster approach, phase transition is modelled by the formation of unstable molec- 
ular aggregates (clusters) whose density is close to that of the new stable phase. This 
argument is the cornerstone of both the equilibrium and the kinetic view of nucleation; both 
of them consider that the aggregation of molecules requires energy, which is known as work 
for cluster formation W{N), and such magnitude can be calculated in terms of the Gibbs 
Free Energy (e.g. chapter 3 of cite 6). The mathematical expression for the work of cluster 
formation is very simple in the case of spherical clusters since the surface tension can be ap- 
proximated by the one known between the phases involved in the process (e.g. vapor-liquid, 
liquid-solid) and does not depend on cluster size. Such approximation is also called capillary 
approximation. In the Classical Nucleation Theory (CNT) the work of cluster formation 
plays a prominent role because the equilibrium cluster size distribution is of Boltzmann 
type, i.e. it is exponentially proportional to the work of formation, and satisfies the law of 
mass action. One of the main problems of this theory is that the pre-exponential factor de- 
pends on the concentration of potential nucleation sites, an obscure parameter which cannot 
be determined precisely. Although CNT is easily comprehensible, in many cases it does not 
accurately predict the nucleation rat^. In the case of Kinetic Nucleation Theory (KNT) 
the pre-exponential factor is derived in terms of the monomer attachment and detachment 
rate^. Although the expression for the collision frequency is readily obtainable, the same is 
not true with the monomer desorption rate. The work of cluster formation enters in KNT 
via the detailed-balance condition imposed to eliminate the dependency on this unknown 
quantity. Considering this assumption one obtains a Fokker-Planck equation as the time- 
evolution equation of the cluster size concentration, better known as the Zeldovich-Frenkel 
equation (e.g. chapter 9 of citCiGJ. Hence, the equilibrium cluster size distribution can be 
re-written in terms of the work of formation and the monomer attachment rate yielding the 
kinetic counterpart of the classical equation. In this framework, the estimated values for 
nucleation rate are also far from those obtained experimentally. Theses deviations led to 
the development of new theories such as the Density Functional Approach. 
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Under DFA, the state of the system is described by the number density of molecules 
as a function of the space coordinates, p(r). No explicit boundary between the phases of 
higher and lower density is assumed and, hence, it cannot accurately be established whether 
a molecule belongs to the old or the new phase. This is also known as the diffuse interface 
theory of nucleation. In this regard, the work of cluster formation is now expressed in terms 
of the Landau potential which can be calculated by using, for example, the square-gradient 
approximationP and used as a replacement of W{N) in CNT and KNT master equations. 

Until recently, the above descriptions of nucleation have been the main tools to predict 
experimental results. Notwithstanding, the underlying mathematical treatment in these the- 
ories is completely deterministic, though they describe a random process. For this reason a 
new formulation for the nucleation theory has been developed within the statistical mechan- 
ics frameworl j^° * ^^ l However, this formalism cannot be used in the case of non-stationary 
conditions and, furthermore, it does not fit to the observed values of (nucleation) induction 
tim(P. The main motivation of this paper is to overcome both limitations, i.e. the deter- 
ministic nature of previous models and the lack of solutions for non-stationary conditions. 

A. Outline of the model 

In the middle of the twentieth century, Hendrik Anthony Kramers described thermally 
activated escape from a metastable state as Brownian motion of a fictitious particle along a 
reaction coordinate, which covers all relevant degrees of freedom that define the system state, 
in a (static) field of force^^. Such processes are characterized by the presence of two stable 
states separated by a potential barrier and can be described in two different but equivalent 
forms. The first is using Stochastic Differential Equations (SDEs) to characterize the time 
evolution of the reaction coordinate. The conjectured SDK is a Langevin-like equatiorP^l 
which contains both a deterministic and a random part. The second is using the time- 
evolution equation of the probability density function (PDF) associated with the reaction 
coordinate. Indeed, this expression can be written as a Fokker-Planck equation^^ (FPE). 
Thus, as we will demonstrate later, using the FPE one can derive an analytical expression for 
the nucleation (or escape) rate in the form proposed by Arrhenius. Moreover, considering 
the nucleation process as a rare event one can obtain the PDF for the induction time. 
Hence, the proposed model not only allows to reduce the problem from a N-dimensional to 
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an unidimensional description but also provides analytical expressions for magnitudes such 
as nucleation rate under both stationary and non-stationary conditions, which demonstrates 
the abilities of the model. 

The aim of this work is to apply Kramers' reasoning to the problem of nucleation in order 
to predict the induction time statistics for both stationary and non-stationary conditions, 
the latter assuming the semiadiabatic limitP'^. We propose cluster radius, R, as a reaction 
coordinate for crystal nucleation and, hence, a stochastic differential equation is postulated 
for modeling time evolution of the cluster size. In this context, analytical expressions for 
statistics of induction time can be deduced by employing the theory of stochastic processes. 
This contribution allows us to better understand the randomness of experimental results 
and the deviation from the theoretical (classical) nucleation rate. 



II. EQUATION OF MOTION FOR THE CLUSTER RADIUS 

A. Theoretical background 

According to the CNT, the work for homogeneous formation of a spherical cluster, i.e. 
the free energy needed to form a spherical cluster, of radius R can be expressed a^ 

AG{R) = -%^R^ + AnaR', (1) 

where A/i is the difference between the chemical potential of the liquid and that of the solid 
phase, Vo is the molecular volume and a is the surface tension between the old and the new 
phase. At small R, the second term of equation ([T]), which is related to the formation of the 
surface, prevails while the first term plays a decisive role at large values of R. The maximum 
of AG defines the critical size, 

with 

(3) 

As widely known, equation ([T| characterizes the equilibrium state of the system whereas 
the cluster radius describes the dynamic of the nucleation process. The latter magnitude 
is generally strongly coupled with the environmental degrees of freedom and due to such a 
coupling its dynamic is not deterministic but stochastic. Therefore, one can assume that 




FIG. 1. Free energy barrier AG with a metastable state in Xji = 0. Escape occurs via the forward 
rate . Red balls represent spherical clusters of the new phase with radius R and black spots 
represent different values of the reaction coordinate X/j which follows a Brownian motion immersed 
in an external field of force derived from AG. The absorbing wall has been denoted as XZ- 



R plays the role of a reaction coordinate in the Kramers' theorjffSHSI, Accordingly the 
variable cluster radius, hereafter denoted as Xji, can be understood as a stochastic process, 
i.e. Xfi{t) draws a random trajectory in the reduced phase spacd^ of the system (figure 
[T|. Considering these assumptions, the nucleation process is characterized by a stochastic 
variable which can be interpreted as the instantaneous position of a fictitious Brownian 
particle. Inspired by the Kramers' theory one postulates that the equation of motion for 
such a particle will be given by a Langevin-type equation. In fact, nucleation is nothing but 
a thermally activated escape process from a potential barrier. Therefore, in order to escape 
from the well A = {Xji < X^}, the fictitious "random walker" must acquire energy to 
overcome the energy barrier and subsequently it must again lose energy to become trapped 
by the attractor B = {Xfj > X^}, i.e. a molecular cluster will fluctuate in size until escape 
from the well A to the attractor B after which the cluster will grow in a deterministic 
manner. 
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B. The model 



The stochastic dynamic of the cluster radius, Xr, can be phenomenologically postulated 
by means of an overdamped Langevin equation, 

r^sMt) = + ./2^m. (4) 

where rj£ is the friction coefficient associated with the reduced phase space (so-called 
reaction-coordinate viscosity) and ^(t) denotes zero- mean, delta-correlated Gaussian white 
noise (GWN), 

m) = 0, 

mm) = Kt-t'). (5) 

It is worth to note that the friction coefficient rj^ is an abstraction that provides information 
about the viscosity of the phase space, £. Hence, there is no trivial expression for rj^ in terms 
of the friction acting on individual molecules in the real space. Nonetheless this magnitude 
can be estimated, together with the mass associated with the reaction coordinate mx-, by 
using molecular dynamic^. One can expect that the former will be related to the mean 



monomer attachment frequency for a i?-sized spherical cluster f(Xj\i) since, the first term of 
the right hand side of equation (|4]) informs about the deterministic behaviour of the cluster 
size and the second one only can arise from the unpredictable collisions of monomers with 
the cluster. Indeed, 7]^ should be characterized by the thermodynamic properties of the 
new phase since at a fixed temperature the monomer attachment frequency depends on the 
density of the final state. 

Thus the most remarkable characteristic of equation Q, comparing with the Ginzburg- 
Landau equation (e.g. chapter 4 of cite [22]), is the second term on the right hand side 
which includes not only a random variable but also the temperature of the thermal bath. 
Hence, this term plays the role of a fiuctuating force which comprehends all degrees of 
freedom associated with the environment. Therefore, the equation of motion of the fictitious 
Brownian particle presents two contributions: i) the deterministic force due to the free energy 
potential AG, and the random force ^{t) := \/2rjskBT^{t). 

From the stochastic differential equation (|4]), using the forward Kramers-Moyal expansion, 



it can be demonstrated that the time-evolution equation of the PDF, p{Xji,t), is 

|(A-„i) = 4.M.f«,<) = -3^, (6) 

where the hnear differential operator is the Smoluchowski operator, 

r _ ^ n dAGjXn) ksT d \ 

a special form of the Fokker-Planck operator Cpi^, and j is the probability current 

Equation ^ is also known as the Smoluchowski equation and belongs to the family of 
Fokker-Planck partial differential equations. 



C. Stationary probability density function 

The time evolution of the PDF converges to the stationary solution of equation ^ when 
t — 7- oo. With the aid of this solution it is possible to estimate the probability to find a 
cluster of a given size when the time is much larger than the relaxation time ts- This 
represents the time required to decay within the attractor A 

-1/2 



1 /a^AG 



mx KO'^Xr 



Xr=0 



(9) 



with mx the effective mass corresponding to the fictitious particle. In such e. when 

the time is larger than the time required for the spontaneous decomposition of a subcritical 
cluster, the non-stationary PDF converges to the stationary one 



C I 

where C is the normalization constant (Appendix [A 



1 / AG(Xk), 

<^ } , (10) 



V / I ( (EAG\ 



The identity (10) is completely in accordance with the cluster size distribution of CNT^ 
with ( playing the role of the pre-exponential factor. 

Nevertheless, the stationary PDF is not useful to obtain the nucleation rate since the 
boundary condition > implies that the probability current must be zero at = 0, 
i.e. j(0) = j{XR) = 0, and hence at any value of X/j including X^. Such restriction implies 
that the escape rate must be equal to zero due to the definition of this magnitude, 

k-" := ^, (12) 
n 

with n the stationary probability that the particle has not crossed the boundary > X'^, 

n = [ %st{XR)dXR, (13) 







which is usually approximated by n ~ 1. This fact shows that the stationary distribution is a 
good but unrealistic approximation. For this reason we shall introduce the quasi-stationary 
solution of the Smoluchowski equation. 

D. Quasi-stationary probability density function 

Now we shall assume that AG/ksT > 1 and T is constant. Under these conditions, 
before reaching the stationary state (i.e. ^ t < oo) the Smoluchowski solution remains in 
a quasi-stationary state and, therefore, the current probability is almost time- independent. 
In such a quasi- stationary state, the probability current over the top of the potential is very 
small near to X^ and the time change of the PDF is also very small. Therefore, the small 
value of the probability current is almost independent of Xr, i.e. j{XR,t) ~ j. Hence, 
equation ([S]) can be written as 

knT _Aniv„\ /h^T d 



J 



Ve (jXr 
Considering now the boundary condition p{X^,t) = 0, one readily gets 

fcijJ Jxu 

III. NUCLEATION RATES AND INDUCTION TIME STATISTICS 

We want to calculate the mean time that the system needs to produce a supercritical 
cluster, i.e. when Xr > X^. This magnitude is also called Mean First Passage Time 



(MFPT) and can be easily related to the induction time. In this section we present the 
expressions obtained for the induction time and the nucleation rate under both stationary 
and non-stationary conditions. 



A. Nucleation rate under steady conditions 



From equations (12), (13) and (15) the following expression for the escape rate can be 
derived, 



1 



V£ ^-AG{z)/kBT 



dz , ' e 



dy e 



AG{y)/kBT 



Accordingly, applying the Laplace (or Gaussian steepest-descent) method (e.g 
cite HM the following approximation can be obtained (Appendix [B]) 



(16) 
. page 124 of 



dXl 



The inverse of this escape rate is also known as the the Kramers time or MFPT, 



(17) 



^AG^/kgT 



The nucleation rate and the induction time are given by the expression 



J 



Poo 

Tk 



t 



ind 



Poo — Pc 
1 _ Tk 

JV ~ Ml 



Id ^AGjO) 



R) 



dXl 



^-AG^/kgT 



TTT]' 



cff 



d^AGiX*) 



(18) 



(19) 
(20) 



where poo is the monomer number density, i.e. the equilibrium number of monomer per unit 



volume J^i/V, and rj'^^ = rje/Mi. As can be observed, equation (19) looks like the classical 
expression for the nucleation rate. Indeed, it can be rewritten as. 



J = Po 



7 7 ^~AG*/kBT 



(21) 



10 



with 



1 d^AG{0) \ 



1 



N) 



dxl 



2nkBT 

being the number of molecules inside the i?-sized cluster, the curvature of the energy land- 
scape at the basin and the Zeldovich's factor (typically® 0.01 < < 1), respectively. 
Therefore, the term ^'^^^'^ ^^^^ time scale of the phase transition. In fact, 

following the classical reasoning (that the characteristic time for nucleation is determined 
by the attachment rate of monomers to the critical cluster) one could postulate that 

^/2k^ fdXN 



so that 



in) - fin) (22) 



% AT7=^ l^i^R) (-23) 



Thus, substituting equation (23) into (21) one gets 



J~Po,ZoZz,/(X;^)e-^^*/'=-^. (24) 



Equation (24) includes the same ingredients as the classical approximation but contains an 
additional correction factor Zq with respect to the classical expression. This new factor is 
of the order of Z^,. Nonetheless, in this work we propose a similar expression for t]^ but 



considering the average value of the collision rate that an individual cluster feels, i.e. /(Xjv). 
Although this magnitude is completely unknown we can assert that it should be lower than 
f{X^). Therefore, our proposed expression for the nucleation rate would be as, 

J ~ Poo ^0 ip^e-^^*/'^^, (25) 



which should predict lower values than the classical one since, /(Xjv) < f{X^). Nonethe- 



less, under this assumption the nucleation rate depends on an unknown magnitude, /(X^v). 
This is the reason by which t]£ must be fitted to experimental values for t-md- Actually, one 



could use the fitted value of rjs in order to estimate the effective value of fiX^) and, then, 
better understand the kinetics of the phase transition. With the aid of such a value infor- 
mation on the underlying mechanism of cluster formation could be obtained by comparing 
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with theoretical values obtained using the different expressions for the monomer attachment 
frequencies (chapter 10 of citeE]). 



B. Induction time statistics under steady conditions. 

Once an estimation for the induction time has been obtained, it seems interesting to derive 
an analytical equation for the induction time statistics which would be helpful for a better 
description of the random nature of such magnitude. Thus, if nucleation is considered as a 
homogeneous Poisson process characterized by the escape rate, A;+ = I/t/^, the induction 
time can be considered as a Gamma-distributed random variable with a PDF given by 

^,(t) = ^exp(-^tl, (26) 



so that. 



P{tind <t) = [ ^i^)ds = 1 - exp {-^i} • (27) 
The latter equation is also called the Kramers Law of the escape time statistics. According 



to equation (18), if we know the viscosity of the reduced phase space rjs then the statistics 



of the nucleation process can be estimated from equation (27). 



Although equations (18), (26) and (27) are good approximations, one can compare the 
accuracy of these expressions by numerical integration. As in the case of ordinary differen- 
tial equations, there exists a huge number of techniques ( Stochastic Integrators , SINT) for 
integrating a SDE as the Euler-Maruyama methocP^. In such a case, the upper limit of Xji 
(i.e. the absorbing wall X^) must be fixed at a value greater than and obeying 

\AG{X^) - AG'(X*)| > 2kBT. (28) 



C. Time-dependent nucleation rate and induction time statistics under the 
semiadiabatic limit 

Let us consider the case when the potential barrier AG changes over time but this change 
is slow compared to the relaxation time of the system r^. According to the reasoning 
of Talkner and co-worker^ the system reaches a quasi-stationary state instantaneously 
(semiadiabatic approximation) and therefore an analogous deduction to the developed for 
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equation (27) can be made for unsteady conditions of AG 



AG*{t)/kBT 



or equivalently, 



k^t) 



/ e2AG(0,t) 


02AG(X*(t),t) 




9X1 


e2AG(0,t) 


a2AG(X*(t),t) 


9X1 


9X1 



(29) 



'AG*{t)/kBT 



(30) 



also called instantaneous escape rate. Proceeding in a similar manner to the previous sec- 
tion, one calculates the expression for the instantaneous nucleation rate and, hence, for the 
instantaneous induction time 



J{t)=p^{t)- 
1 



Vs 



^ind(^) 



J{t)V Mi{t)' 



(31) 
(32) 



Equation (29) constitutes a good approximation only when the topology of the potential 
barrier AG(Xij,t) does not change, i.e. the character of its local maxima and minima 
must be invariant. In fact, Talkner and co-worker^ derived a more accurate expression 



which contains expression (29) and a second-order correction term. They called the latter 
as geometric correction term because it is related to the geometric change of the barrier 
shape. We will not consider this second-order correction in this work. 

Following a similar procedure to that of the previous section, the nucleation process can 
be understood now as a non-homogeneous Poisson proces^ and, hence, the first-passage 
time statistics is governed by the equation below 



with 



and assuming 



git) = P^e^v{-RK{t)} 

RK{t) 



Po 



-ds. 



Rk(oo) 



'rK[S 
Poo(s) 



ds 



oo. 



Accordingly, the distribution function of the induction time is given by 



P{Und<t)= [ f?(s)ds = l-e-^W. 
Jo 



(33) 

(34) 
(35) 

(36) 
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The accuracy of our theoretical predictions is strongly related to the assumptions made 



above. Consequently, the results predicted with the aid of equations (29) and (36) can 
be improved using a SINT. In fact, the stochastic integration is indispensable in order to 
obtain more realistic predictions when the hypotheses of the semiadiabatic limit are not 
fulfilled. Indeed, SINTs are the only tool to predict escape rates and, then, induction times 
under strongly unsteady conditions. As mentioned in the previous section, the instantaneous 
position of the absorbing wall (t) must be fixed at a value greater than X^(t) and obeying 
the following relation 

|AG'(X^,t) - AG{X*^,t)\ > 2kBT. (37) 



IV. THE MODEL AT WORK 

The Brownian model proposed in sections 2 and 3 has been tested by fitting the viscosity 



parameter through the expression for the induction time (20) to experimental results of such 
a magnitude measured for hydrated calcium sulphate CaS04-2Il20 (gypsum) in a volume 
of 200 fiL"^^. This kind of applications represents one of the main contributions of this paper 
because neither classical nor non-classical nucleation theories allow to follow a nucleation 
event in such a simple manner. As a summary, one must fit rjs to experimental data using 
equation (20) and apply both SDE ^ and equations (19)-(27) and (31)-(36) to predict both 
nucleation rates and induction time statistics under both steady and unsteady conditions. 



Once we have the estimation of r]£, one can apply equation (27) in order to calculate the 
induction time statistics and verify the ability to predict P{tind < t). Moreover, with the 
aim of considering a finite number of assays the time evolution of the reaction coordinate X^ 
was simulated using the Euler-Mauryama method. Furthermore, using the fitted value of 
rj^ one could simulate more complicated experimental conditions where classical expressions 
cannot be used. 

Twelve assays measuring induction time were considered at a fixed supersaturation and 
temperature, e.g. S = C/Ce = 1.89 at T = 328K, to estimate r]£, resulting in a value r]'^^ ~ 



4.511^. When substituting such a value into equation (18) one can calculate a theoretical 



curve of predicted induction times as a function of supersaturation. The predicted curves 
for induction times are in very good agreement with the measured values (figure [2]). 

The model also provides the necessary tools for calculating the induction time statistics by 
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FIG. 2. Experimental and theoretical values of i^rf as a function of supersaturation C/Ce at 
T = 328 if (left panel) and T = 331 if (right panel). Dashed lines correspond to induction 



times predicted by using equation (18) and the estimated friction coefficient •q'^^ ~ 4.51. Each 
experimental value represent an average of twelve assays. 



considering the Kramers law (eqs. ( 27 ) and ( 36 ) ) . The same set of twelve experimental values 



for induction time were used in order to calculate the experimental cumulative distribution 
function (blue triangles in figure^. Using the estimated value of r]^^ into equation (27), 



Pitmd < t) (solid orange hne in figure [3]) was obtained. As can be observed, the theoretical 
prediction does not fit as expected to the experimental statistics. Such a disagreement 
between the predicted and the observed curves is not due to the assumptions made to 
integrate the escape rate but to the fact that in experiments we truncated the statistics to the 
longest time observed, ^^£0 i-^- nucleation events were recorded after an arbitrary time 
^max corresponding to the duration of the experiment. So the experimental tind values were 
biassed towards small values. Using a SINT (Euler-Maruyama) to simulate with equation 
^ a finite number of assays with an upper limit equal to the experimental observation time 
(i.e. if a simulation exceeds the upper limit tmax, then such assay is not considered and 
another one starts) the simulated results (green squares in figure |3| are much closer to the 
experimental curve. 
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FIG. 3. The cumulative distribution P{tind < t) as a function of time: i) experimental data (blue 
triangles) and ii) simulated results (green squares). Solid orange line represents the cumulative 



distribution function given by equation (27). 
V. CONCLUSIONS 



In this work a Stochastic Differential Equation was presented as the equation of motion of 
the classical order parameter in the Nucleation Theory, i.e. the cluster radius. This SDE for 
modelling the time evolution of the radius of an individual cluster was used to obtain theo- 
retical equations to predict (nucleation) induction times and its statistics. These equations 
are applicable to the often used crystallization setups in which supersaturation changes over 
time (e.g. due to changes in bulk concentration, temperature, pH, etc.^^) and can be used to 
identify the underlying mechanism of cluster formation by fitting the measured nucleation 



rates to equations (29)-(29) and (25), as well as to predict induction times (eqs. (20) and 



(32)) and induction time distributions (eqs. (26)-(27) and (33)-(36)). We present here not 



only an application of the model but also a method to obtain theoretical and simulated 
predictions of both induction times (and hence of nucleation rates) and cumulative distribu- 
tion functions at different concentrations and temperatures. The first test of the stochastic 
model against experimental data reveals its potential ability for calculating nucleation rates 
and induction time statistics. Both the analytical and the numerical results predicted by 
the Brownian model seems to be in good agreement with the experimental data. Hence, 
the interpretation of the nucleation stage as an escape process could be an optimal tool to 
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study deeper problems of the first-order phase transitions. However the theoretical approx- 
imations usually offers worse results than the simulations of the SDE by using stochastic 
integrators due to the experimental truncation of the statistics to the longest time observed. 
Additional work is in progress to deduce analytical expressions for nucleation rate, induction 
times and induction time distributions when neither stationary nor semiadiabatic limit can 
be considered. 
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APPENDIX 

Appendix A: Approximation of the normalization constant 

Let us consider that i^G* /ksT is large and that ksT is very small. Therefore pst becomes 
very small for values of Xr appreciably different from Xr = 0. In this case AG can be 
expanded according to Taylor's theorem as 




(Al) 
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and thus one gets the following approximation 



exp 







AG(0) + ^AG"iO)X^ 



-AG{0)/kBT 



-AG{0)/kBT 



B-' 

oo 



dX 



^AG"(0)X 



AG"(0)X2|rfX 

X=oo 



2 \ v^2%T 



x=o 



with 



^^g-AG(0)ABT 

2 VAG^"' 



AG'"(X) = 



g'AG(X) 
9X2 



denoting the second derivative with respect the reaction coordinate X. 



Appendix B: Integration of the escape rate equation 



(A2) 



Whereas the main contribution to the first integral in equation ( 16 ) stems from the region 



around = 0, i.e. close to the minimum of the barrier, the main contribution to the second 
integral stems from the region around X^. Therefore, considering Taylor's expansion of AG 
around its minimum and maximum, 



x{Xn-0r + O{Xl), 



AG{Xr) ^ AG(X* 



1 [d'^AG{X*R] 



2! 



X 



X (XH-X*)2 + 0((Xfi-X 



R) 



(Bl) 



(B2) 



and substituting equations (Bl) and (B2) into equation (16), one obtains the escape rate 
equation. 



k+ 2 uquj* 



erf 



co%Xn-X} 



XI 



(B3) 
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with 



Wo 



In the case that > 1 and (X^ — X^) > 1 one can consider 



erf 



Xr=0 



as a good approximation. Finally, one readily obtains the desired expression for the escape 
rate 



(B4) 



Nevertheless, the accuracy of the result given by equation (B4) depends on the goodness of 
the approximations assumed above. 
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